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m I Abstract. We describe systematic changes in the centroid frequencies 

' and the sphtting coefficients as found using data from MDI on board 

CN i SOHO, covering cycle 23. The data allow us to construct a seismic map 

of the evolving solar activity - covering all latitudes. At lower latitudes, 
the temporal evolution closely tracks that of butterfly diagram. The ad- 
' ditional information from higher latitudes in the map is of a significant 

! activity in the polar region, peaking at activity minimum in 1996. The 

' most plausible source of solar oscillation frequency changes over the solar 

cycle is the evolution of the radial component of the small-scale magnetic 
field. The amplitude of the required mean field changes is ~ 100 G at the 
photosphere, and increasing going inward. 
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2 ■ !• Phenomenology of frequency variations 

, Using satellite intensity data from cycle 21, Woodard & Noyes (1985) first noted 

that p-mode frequencies increase with increasing solar activity. This discovery 
was confirmed by a number of investigations made during cycle 22. In par- 
^ ' ticular, it was then established that the most significant changes occur in the 

. antisymmetric part of the spectrum of oscillation multiplets - that is, the part 

that refiects the asphericity of the sun (Kuhn, 1988; and Libbrecht &; Woodard, 
1990). The dominant part of the multiplet structure is symmetric and arises 
from advection by rotation. 

Oscillation data, from the on-going spatially resolved observations, are rep- 
resented as centroid frequencies, and splitting coefficients, a2k,en, in the 
following expression for the frequencies of the individual modes 



l^inm — ^en + '^ 0-k/n'Pk,t- (1) 
fc=l 



The quantities V are orthogonal polynomials of m defined for 2k < i (see Ritz- 
woller and Lavely 1991 and Schou, et al. 1994). 

We use MDI data, which contains coefficients up to A: = 18 for about 2000 
p- and f-modes with i < 300. There are 24 sets of data corresponding to 72 day 
long measurements done between May 1996, when the sun was at its activity 
minimum, to June 2001, when the sun was in it high activity phase. 
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The relative changes of solar frequencies are of order 10~^, which does is 
comparable to the individual measurement errors. Significant rates of change are 
obtained by binning the data or by assuming, as Libbrecht and Woodard (1990) 
did, that the changes scale inversely to mode inertia, which what is expected if 
the activity related changes acts only near the surface. Following this idea, we 
write the frequency changes in the form 

^e,n - {i^e,n)niin = 7^, (2) 

where the subscript "min" refers to measurements made at solar minimum. We 
assume the same scaling for the variable part of the a coefficients of even orders, 
and thus, we write 

0-2k,e,n = 0'2k,e,n;y:ot + Ck,ej — , (3) 

The numerical factor Ck/ (see Dziembowski, et al, 1999, for the explicit expres- 
sion and its justification) was introduced to make each 7^ an unbiased probe of 
the distortion of the sun. Here, we removed the constant contribution from 
the centrifugal force, which is the only non-negligible effect of rotation. 

The constant values of 7^ inferred from the p-mode data show systematic 
evolution as solar activity varies. The coefficients, 7^, up to /c = 9 are in excellent 
correlation with the corresponding coefficients from the Legendre polynomial 
expansion of the Ca II K line intensities, which are regarded as a good proxy 
for the magnetic field (Dziembowski et al. 2000). 

For p-modes, we expect that the 7's do not depend on £, if their source is lo- 
calized near the surface. A priori, we might expect a significant frequency depen- 
dence, but it is not very strong. This tells us something the precise localization 
of the source. For the f-modes, where we have the approximate proportionality, 
ly oc a component of Au was determined, which is also oc y/I and grows 
with increasing activity. This component may be interpreted as evidence for a 
contraction of the sun's outer layers as activity rises. Such contraction would 
take place if the increase of the field were dominated by its radial component 
(Dziembowski. Goode & Schou, 2001) 

2. Seismic map of the sun's activity 

Having determined the 7^ coefficients as functions of time, we can construct a 
seismic map of varying solar activity. To do this, we determine the quantity, 

j{e,t) = Y,ikit)P2kicose). (4) 

fe=0 

In Figure 1, we show the changes in 7 averaged into bins 5°-wide in latitude. 
At lower latitudes, we recognize features that are well-known from the butterfly 
diagram: activity appearing first at 35 — 45°, and gradually moving toward the 
equator. However, there is also unexpectedly large activity in the polar regions. 
This decreases with the rise of activity at the lower latitudes. This result is not 
an entirely surprising one, because it has been known for some time that the 
polar field flips at activity maximum. Furthermore, Moreno-Insertis &; Solanki 
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Figure 1. Zonal averaged values of 7 in //Hz. The ranges of the 
latitude are given in the legend, and the overall zero-point is arbitrary. 



(1999) have already found that their models for low i mode frequency behavior, 
during the activity cycle requires, fields at high-latitude. 

3. The cause of the frequency changes 

There is ample evidence that the behavior of the 7s in some way reflects changes 
in the magnetic field in the outermost layers of the sun. We would like is to infer 
more about the field: its intensity and its behavior as a function of the depth 
below the photosphere. For this it is essential to consider the role of changes 
in the temperature and turbulent velocity induced by the magnetic field. A 
significant role of the latter effect may be eliminated (Dziembowski and Goode, 
in preparation). 

Following Goldrcich ct al. (1991), we consider a small-scale, random mag- 
netic field as the primary cause of the frequency changes in which we allow the 
vertical (r) component to be statistically different from the two horizontal {6 and 
(f)) ones. In this, the averaged values Bj are treated as functions of depth, and 
slowly varying functions of the co-latitude, The latter dependence is represented 
in the form of a truncated Legendre polynomial series, 

BiBj = 6ij^[SjrMr,k{^) + ^MH,ki^je + Sj^]P2k{cos9), (5) 

A:=0 

where we included only seismically detectable (symmetric about equator) terms. 

Each of the ^-components gives rise to a distortion of the sun's struc- 
ture. For k > 0, the hydrostatic equilibrium condition suffices to determine the 
distortion of all the thermodynamical parameters. 
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In the evaluation of the 7s, which in general must be treated as functions 
of both v and £, we use the variational principle for stellar oscillations in which 
we treat the effects of the magnetic field as a small perturbation. This principle, 
with use of hydrostatic equilibrium and equations for adiabatic oscillations, leads 
to 

lk = J (l^k,T^ + IC^,rMr,k + }Cj^^HMH,k^ dD, (6) 

where D is depth. All the kernels, /C, may be explicitly expressed in terms of 
parameters of the standard solar model and the radial eigenfunctions of its p- 
modes. If the magnetic perturbation is significant only in the layers well-above 
the lower turning points of all p-modes considered, then the kernels are both k 
and £ independent. 

Goldreich et al.(1991) considered only changes in centroid frequencies. They 
pointed out that to explain the frequency increase between 1986 and 1988 a 1% 
decrease of photospheric temperature is needed {IC^^t is always < 0). Regarding 
this requirement as being incompatible with observations, they adopted the 
changing magnetic field as the sole cause of the frequency increase; they found 
that the field increase must be about 250 G at the photosphere, and steadily 
growing to about 1 kC at a depth of 10 Mm. Their numbers refer to the case 
of a statistically isotropic field {MH,k = '^■^r,k)- A much more modest field 
increase (< 100 G at the photosphere) would result for an inwardly growing, 
pure radial field. A similar result was obtained by us (Dziembowski, Goode k, 
Schou, 2001) from our analysis of the MDI data. 

For A: > 0, the term involving temperature change may be eliminated. In the 
outer layers, the signs of the modified kernel are > and /C^^ < 0. Thus, 
the most economical requirement is a pure radial field. With the assumption 
that all the frequency changes are due to an increase of the radial component of 
the small-scale magnetic fields, the plots in Fig. 1 may be scaled from fiHz to 
Gauss using the factor ^ (190G)V1a*Hz. 
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